Validity of Anderson's theorem for s-wave superconductors. 
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We investigate validity of Anderson's theorem (AT) for disordered s-wave superconductors in a 
negative U Hubbard model with random on-site energies, e». The superconducting critical temper- 
ature, T c , is calculated in the coherent potential approximation (CPA) as a function of impurity 
concentration, c, the random potentials for different band filling. In contradiction to Anderson's 
theorem, we found that T c is dramatically sensitive with respect to c and Si. Our results shows 
that for low impurity concentrations and weak on-site energies, ei, the AT is valid, while in the 
strong scattering limit even for low impurity concentration, T c is very small with respect to the 
clean system and by increasing c it is completely suppressed, hence AT is violated in this regime. 



I. INTRODUCTION 

The insensitivity of the metallic conventional su- 
perconductors to nonmagnetic impurity doping is well 
known and was be explained in terms of Abrikosov and 
Gorkov (AG) theorjii and Anderson's theorem (AT)—. 
According to this theorem the superconductor critical 
temperature, T c , and the superconductor density of 
states are not affected by the nonmagnetic impurity 
scattering. In contrast to this behavior in conventional 
metallic superconductors, in the high temperature 
superconductors (HTC) doping of a small amount 
of nonmagnetic impurities like Zn destroys super- 
conductivity, as reported for Y Ba2(Cu\- x Zn x )j,0-, ? ' , 
YBa 2 (Cui- x Zn x ) i O g 2£, La\.mSros5Cu\- x Zn x O£, 
Bi2Sr2Ca(Cui- x Zn x )20gL. Physical evidence such 
as T 2 dependence of penetration depth in disordered 
HTC 8,9 and STM measurements on single Zn doped on 
Bi2Sr2Cu20g,+s^, shows that non-magnetic impurities 
behave as strong scattering centers. While this can 
be understood as a consequence of d-wave pairing in 
HTC, this kind of strong-scattering impurity role has 
raises the question for strong scattering effects in s-wave 
superconductors? In particular, do the order parameter, 
density of states and especially T c change with impurity 
doping. In other words, is the Anderson theorem valid? 
Ghosal et alii calculated the local order parameter and 
the density of states of a disordered s-wave supercon- 
ductor with a uniform impurity distribution, [— V, V]. 
They found that in the strong scattering limit, V — * oo, 
for some regions of space the local order parameters 
are zero while in some others they are finite. A binary 
alloy s-wave superconductor has been investigated by 
one of us—, and it has been found that in the strong 
scattering limit the average order parameter goes to 
zero and a superconductor to insulator phase transition 
is take place, so the AT is violated in this case. Also 
the existence of an impurity bound state in a disordered 
s-wave superconductor energy gap is discussed by 
Balatsky et alm&. 

In this paper we investigate the validity of AT for s- 
wave superconductors in detail, by calculating T c as a 
function of impurity concentration and electron band fill- 



ing. This is in contrast to previous work challenging AT 
in s-wave superconductorsi^ii^ which were calculations of 
the order parameter and density of states, not directly of 
T c . 

We found that at low impurity concentrations and in 
the weak scattering (gj (( band width) regime, the devi- 
ation of T c with respect to the clean system critical tem- 
perature T c q is negligible, and hence AT is satisfied. But 
in low impurity concentrations at the strong scattering 
limit, although superconductivity is preserved, the T c is 
much less than T c q, so the AT is not obeyed. By increas- 
ing c in the strong scattering limit T c is completely sup- 
pressed hence a superconductor-insulator quantum phase 
transition take placed^. Furthermore, we also calculated 
T c at a fixed averaged band filling for different impurity 
on-site energies and impurity concentration. Our results 
show that T c is decreased by increasing impurity concen- 
tration and eventually goes to zero in the strong scatter- 
ing regime. 

II. MODEL AND FORMALISM. 

We start our investigations with the following at- 
tractive Hubbard model on a two dimensional square 
latticed, 

H = - *ijC- .Cj -+y^ Uihigni-g+y^^Si-iJ^hia. (1) 

ij ia ia 

where tij are the hopping integrals, Ui is the local attrac- 
tive interaction between electrons of opposite spin, and 
H is the chemical potential. Finally, Ei is the random on- 
site energy which takes values of with probability c for 
impurity sites and S with probability 1 — c for the host 
sites. 

In the mean field approximation, for any given config- 
uration of impurity sites the solution of Eq.^is given by 
the Bogoliubov de Gennes equation 
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where H a = (gj - Ui)8 a - t$ and fa = fj, - \ll l n l is 
the on-site renormalized chemical potential due to the 
usual Hartree decoupling. The local order parameter Aj, 
assumed to be real, is defined by, 

/oo 
dE f{E) ImG 12 (i,i;E) (3) 
- OO 

and the local band filling n t ; is, 

/oo 
dE f(E) ImG n (i,i;E) (4) 
-oo 



where f{E) 



—ip is the Fermi function. Eq[2] 



can be written as, 
G(i, j; £) = G°(i, j; BJ+^G^, I; £)V;G(J, j; E) 



where the random potential matrix V; is 
ei + Wini A/ 



V, = 



A ; * -e, - i£7,n z 



(5) 



(6) 



and the clean normal system Green function G (i,j;E) 
is given by, 



£-e k 
£ + e k 



. (7) 



„k.r; 



is the bulk band struc- 



where e k = -fj, - jf Y,ijt\_ 
ture. In our numerical calculations hopping to the first 
and second nearest neighbors is considered and we ne- 
glected other hopping terms, so e k = — u — 2t(cos(k x a) + 
cos(k y a))+t cos{k x a)cos{k y a) . The interaction potential 
is chosen to be [7, = 2t. 

The Dyson equation for the average Green function 
,G(i,j;E) , corresponding to EqJSJis, 

G(i,j;E) = G°(i,j;E) 

+ J2 G °(i>hE)V(l,l' ,E)G(l ,j;E), (8) 
u' 

where the self energy, E(£, I ,E), is defined by, 

(V«G(i, j; E)) = £ S(Z, l,E)G(l',j; E). (9) 



The Fourier transformation of Eq[S]leads to the following 
relation for the average Green function G(i,j;E), 

G(i,j;£) = ^E k e* k - r « 

E-ek-Eutk,^) -E 12 (k,£) 
-E 21 (k,£) £ + e k -E 22 (k,£:) 



(10) 

Ea l3llUI should be solved to determine the averaged 
Green function, G(i,j;E). In general there is no 



analytical solution for such systems, so it should be 
solved approximately. Many approximations such as self- 
consistent Born approximation or Abrikosov-Gorkov the- 
ory where is valid for weak scattering (<5 (( band width) 
and all impurity concentrations^, self-consistent T- 
matrix approximation where is limited to the case of low 
impurity concentrations 14 and the coherent potential ap- 
proximation (CPA) where is valid for all impurity con- 
centrations and impurity strengths 15 , S, and also recover 
both self-consistent Boron approximation (SCBA) and 
self-consistent T-matrix approximation (SCTMA) 1 ^ 1 ^, 
have been applied. To compare Abrikosov-Gorkov the- 
ory (weak scattering limit) and strong scattering case 
(5 )} band width) we use CPA where in the weak im- 
purity strengths (S {{ band width) is equal to the self- 
consistent Boron approximation of Abrikosov-Gorkov 
theory. 



III. CALCULATION OF T c IN THE CPA 
FORMALISM. 

In the CPA formalism inter-site correlations are ne- 
glected and each lattice site is replaced with a site of an 
effective medium except one which is called the impurity 
site and is denoted by i. So the self energy is local and 
takes the same value for all sites, £(/, m, E) = H(E)Si m , 
hence Easl^landHUIare reduce toi^, 



(V i G im v(i,i;E)} = -£(E)G(i,i;E) 



(11) 



and 



G(/,m;£;) = l^ e k - 



E-e k -X n {E) -E 12 (£) 
-E 21 (£) £ + e k -E 22 (£;) 



(12) 



respectively, where the impurity Green function, 
G lmp (i,i, E), is related to the averaged Green function 
G(i, i, E) as±i, 

G imp {i,i;E) = G{i,i;E) 

+ G{i,i;E){V t -V{E))G m v{i,i-E). 

(13) 

The new average Green function G(i,i;E) is given by 
taking an average over all possible impurity configura- 
tions, 



G(i,i;E) = (G im "(i,i;E)) 



(14) 



Easllll- IT41 should be solved self consistently to provide 
G(i,i;E). We continue our discussions by derivation of 
a set of equations to obtain the critical temperature in 
the CPA formalism. 

At T c both the local and average order parameter Af 
and A go to zero, and the spatial deviation of the local 
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order parameter, <5A;, with respect to the averaged order 
parameter, A, is negligible, so 

Ai a A. (15) 

By inserting Ea 1 151 in Ea II II we found that, 

((e i --U i n i )G n (i,i;E)) + AG 21 (i,i;E) = 

^ 11 {E)G{i,i;E) + X 12 (E)G 21 (i,i;E).(l6) 

Since at T c the system is normal, the first terms in the 
left hand side and right hand side of Eq^Jare equal and 
cancel each other. Therefore we find, 



E 12 (£) = A. 



(17) 



By inserting EqElinto the gap equation, EqOl for the 
average order parameter, A, at T c we obtain, 



/oo 
dEf(E) x 
-OO 

Im ( f de 

J — OO 



(E-e- Z n (E)(E + e - E 22 (E)) 



) 

(18) 



where No(e) is the clean system density of states defined 
by, 



ATo(e) = ^5>(e-e k ). 



(19) 



Eqs^l an( l [HI Eas llllTil should be solved self consis- 
tently to calculate En(-E) and E 22 (-E). Then these self 
energies should be inserted in to Eq^|to obtain the crit- 
ical temperature, T c . 



IV. RESULTS AND DISCUSSION 

We considered several cases, first in the weak scatter- 
ing regime, S — 0.05i where CPA is equal to SCBA ( 
Abrikosov-Gorkov theory), T c is calculated in terms of 
average band filling, n, for different impurity concentra- 
tions. Fig^ shows T c in terms of n at a fixed impurity 
strength, 6 = 0.05i, for different impurity concentrations. 
As one can see, T c is changed only very little with respect 
to the clean system, and hence Anderson's theorem is 
satisfied in weak scattering, as expected. 

Second, T c was calculated in terms of n at fixed impu- 
rity concentrations, c — 0.1, c — 0.15 and c = 0.2, for dif- 
ferent value of 6 varying from weak (SCBA or Abrikosov- 
Gorkov theory) to strong scattering. Fig|21 illustrates our 
result in this case. For all these impurity concentrations 
in the strong scattering regime T c decreased with respect 
to the clean system significantly and for some band fill- 
ings there is no superconductivity. Also we found that 
for high impurity concentrations, c = 0.15 and c = 0.2 
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FIG. 1: T c in terms of average band filling in the weak scatter- 
ing regime (SCBA or Abrikosov-Gorkov theory) for different 
impurity concentrations, c, where S — 0.05t (( bandwidth. 
The deviation of T c from the clean system critical tempera- 
ture, T c o, is small. Hence AT is satisfied in the weak scattering 
and low impurity concentration regime. 



a superconductor-insulator quantum phase transition is 
take placed as previously predicted by usi^ 

Thirdly, T c was calculated for fixed values of the scat- 
tering strengths, S — 2t and <5 — 20t, and different values 
of impurity concentration. We found by increasing impu- 
rity concentration T c is reduced and at a critical impurity 
concentration superconductivity completely disappears. 
FigGD shows the results in this case. 

Finally, at a fixed band filling, n — 0.9, the critical tem- 
perature T c is plotted as a function of impurity concen- 
tration for different value of 8. Fig0|shows that for weak 
scattering S (( band width (SCBA or Abrikosov-Gorkov 
theory), the critical temperature T c is not much influ- 
enced by increasing impurity concentration while in the 
strong scattering regime T c is decreases linearly with im- 
purity concentration and superconductivity is completely 
destroyed. 

In conclusion, we have investigated response of the crit- 
ical temperature T c of an s-wave superconductor to non- 
magnetic impurity doping. Three cases are considered. 
First, in the weak scattering regime 6 (( band width 
(SCBA or Abrikosov-Gorkov theory regime), T c is cal- 
culated for different impurity concentrations. Here we 
found that the T c is not affected by impurity doping so 
the Anderson theorem is satisfied in the case of low impu- 
rity concentrations and weak scattering. Second, at fixed 
impurity concentrations c = 0.1, c = 0.15 and c = 0.2, T c 
is calculated in terms of band filling for different value of 
S. We found that by increasing 5 the critical temperature 
decreased and at a critical S a superconductor-insulator 
quantum phase transition takes place. Third, at fixed 
S — 2t and S = 20t the critical temperature is calculated 
in terms of band filling. Our results in this case show 
that by increasing impurity concentration T c is again re- 
duced, and at a critical concentration a superconductor- 
insulator quantum phase transition again takes place. Fi- 
nally at a fixed band filling n — 0.9, T c is plotted in terms 
of impurity concentration. In this case we found in the 
weak scattering regime (SCBA or Abrikosov-Gorkov the- 
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FIG. 2: T c for c = 0.1, c = 0.15 and c = 0.2 in terms of band 
filling for different 8 varying from weak (SCBA or Abrikosov- 
Gorkov theory) to strong scattering. In all these strong scat- 
tering cases T c is sensitive with respect to c and <5, so AT is 
not valid. For c = 0.1 superconductivity is preserved even 
in 8 = lOOt but for c == 0.15 and c = 0.2 superconductiv- 
ity is suppressed at a critical 8, so a superconductor-insulator 
quantum phase transition is take place. 



ory regime) T c is not much influenced by increasing im- 
purity concentration but at strong scattering regime T c 
is linearly decreased by increasing impurity concentra- 
tion. Therefore, we conclude that Anderson's theorem^, 
only holds in the weak scattering (Born approximation) 
regime, originally considered by Abrikosov and Gorkov 1 . 
In stronger scattering even non-magnetic impurities can 
suppress s-wave superconductivity and lead to greatly re- 
duced T c or even a transition to a non-superconducting 
insulator state. 
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FIG. 3: T c in terms of average band filling at hxed on site 
energies 8 — 2t and 8 = 20/ for different impurity concentra- 
tions. By increasing impurity concentration not only is T c re- 
duced, but also at a critical impurity concentration supercon- 
ductivity is completely suppressed, hence a superconductor- 
insulator quantum phase transition has occurred. 
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FIG. 4: T c in terms of impurity concentration, c, for different 
values of scattering strength, 8, where the average band filling 
is fixed at n — 0.9. For small 8 (SCBA or Abrikosov-Gorkov 
theory regime) T c is not influenced much and AT is valid, 
but in the strong scattering regime T c is reduced linearly by 
increasing impurity concentration until T c eventually becomes 
zero. 
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